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New 13 CO observations towards the Keplerian gaseous 
disk around R Mon 
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Abstract We present new high angular resolution ob- 
servations of the 13 CO 1— >0 rotational line towards the 
HBe star R Mon, obtained with the IRAM Bure In- 
terferometer, along with the previous results in the ob- 
served transitions 12 CO 1— >0 and 2— >1. We have used 
a flat disk model to fit the 12 CO 1— >0 and 2— >1, and 
13 C0 1— >0 emission, in a strip perpendicular to the 
outflow axis. The model assumes standard abundances 
(X( 12 CO) = 8-1CT 5 , X( 13 CO) = 9-10" 7 ), radial poten- 
tial temperature and density laws (T (K) = T Q r~ 9 , 
p (cm -3 ) = p v~ p ), and Local Thermodynamic Equi- 
librium. The 13 CO and 12 CO emission is consistent with 
a flat disk at an inclination angle of 20° in Keplerian ro- 
tation around the star. The gaseous disk is fitted with 
a mass of 0.014 M Q , an outer radius of 1500 AU, a 
temperature of 4500 K at the inner radius (1 AU), and 
values of q=0.62 and p=1.3 for the indexes of the tem- 
perature and density laws. This new 13 CO observations 
allow us to conclude that the disk around R Mon is flat. 
Our result confirms previous works suggesting a pre- 
dominant flat geometry in disks around early Be stars. 
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1 Introduction 

A big theoretical and observational effort has been done 
in recent years for the understanding of the disk occur- 
rence and evolution in Herbig Ae/Be stars. Evidence of 
the existence of dusty and gaseous circumstellar disks 
around some HBe stars exists at optical, near and mid 
IR (see for example lMeeMS et al. 2001 [Wink et al. 200~2~\ 
Millan-Gabet et al. 200 1\ \Acke et al. 2005\ and more 
recently \lsella et al. 2006} . It is generally accepted the 
presence of disks similar to those in T Tauri stars for 
Herbig Ae and late (spectral type later than B7) Be 
stars, usually flared disks with inner rims, but the asso- 
ciation of more massive early Be stars with disks is more 
uncertain, and the direct detection of a circumstellar 
disk around an HBe star at millimeter wavelengths has 
remained elusive until recently. The direct detection of 
a gaseous disk around R Mon also allows us to study 
the geometry and kinematics of a circumstellar disk in 
a B0 star. 

R Mon is a very young B0 star with a T Tauri 
companion separated by 0.69" . Because of its youth, 
R Mon is not directly visible in the optical, but ap- 
pears as a resolved conical reflection nebula in scat- 
tered light. At infrared wavelengths, R Mon appears as 
a point source located 0.06" ±0.02 south from the opti- 
cal peak d Weigelt et al. 2002~\ . Based on HI emission- 
line ratios, Close et al. (1997) measured an extinction 
of Av = 13.1 mag towards the star that they interpreted 
as due to an optically thick disk of R = 100 AU. 



2 Previous results and new observations 



Fuente et al. (2003) detected a dusty disk around R 



Mon. New continuum observations \Fuente et al. 2006\ 
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hereafter Paper I) carried out with the A configura- 
tion of the PdBI allowed us to constrain the size, mass, 
and spectral index of this circumstellar disk. Our re- 
sults showed that the size of the dusty disk is 0.3", 
or 150 AU of radius, in agreement with the value de- 



rived by | Close et al. (1997) Fitting the SED in the mm 
range (See \Fuente et al. 2006\ for details) we obtained 
our best fit with a dusty disk of mass = 0.007M Q 
and a spectral index (3 — 0.3. The low value of (3 sug- 
gested us that grain growth has proceeded to large sizes 
in the disk around R Mon. Typical values of 0.5-1 are 
found in Herbig Ae and T Tauri stars, also associated 
with grain growth. 

Interferometric observations of 12 CO 1— »0 and 2— >1 
transitions were also presented in Paper I. We discussed 
our data by means of the results of a radiative transfer 
model applied to the observed molecular emission. We 
found that the observations were well fitted using a flat 
gaseous disk in Keplerian rotation, with an apparent 
inclination of 20 degrees from edge-on view, and a mass 
of 0.014 Mq. We were also able to constrain the mass 
of R Mon to 8±1 M Q . We found some keys to support 
a flat disk instead of a flared disk geometry. However, 
since the CO emission from this circumstellar disk is 
optically thick, we were not able to discard a flared 
geometry for the disk. 

13 CO 1— »-Q observations were carried out with the 
iramS array at Plateau de Bure between January 2006 
and March 2006, using AB configuration (beam l.l"x0.7" 
at 110.2 GHz). The signal-to-noise ratio (5a detection) 
is not good enough to study the spatial distribution of 
this molecule, but allows us to conclude about the ge- 
ometry of the disk around R Mon. All images in this 
paper are centered at RA 06:39:09.95, DEC +08:44:09.6 
(J2000). 



3 Model results and discussion 

The interferometric CO images (figureQJ; and QJ1) show 
a compact molecular clump towards the star position 
with a size of 2.5" (2000 AU at d = 800 pc) in the 12 CO 
2— >1 line and a more extended and weaker component 
which follows the shape of the nebula. In figure [2] the 
blueshifted and redshifted components are drawn sep- 
arately. The orientation of the major axis of the disk is 
found to be perpendicular to the outflow axis. 

The Position- Velocity (P-V) diagrams (figure[3|) show 
the emission along the major axis of the disk. We have 
modeled the 12 CO 2— >1 and 1— >0 emissions using a 

1 IRAM is supported by INSU/CNRS (France), MPG (Ger- 
many), and IGN (Spain). 
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Fig. 1 (a) and (b) show the interferometric images of the cir- 
cumstellar disk around R Mon in the continuum at 1.3mm (a) 
and 2.7mm (b). Contour levels are 1.5 mjy/beam to 8mjy/beam 
by 1.5mjy/beam in the 1.3mm image, and 0.75 mjy/beam to 8 
mjy/beam by 1.5 mjy/beam in the 3mm image. The images have 
been obtained by convolving the clean components with a 0.3" 
circular beam in the case of the 1.3mm map, and a 0.8" circular 
beam in the case of the 2.7mm image. In (c) and (d) we show 
integrated intensity maps of the 12 CO 2— >1 (c) and 12 CO 1— >0 
(d) rotational lines obtained with the PdBI. The contour levels 
are 0.25, 0.5 to 6 by steps of 0.5 Jy/beam*km/s in the 12 CO 2->l 
image and 0.25, 0.5 to 6 by 0.25 Jy/beam*km/s in the 12 CO l-»0 
image. 



Parameter 


Fitted value 


R Mon mass (Mq) 


8 ± 1 


i (deg) 


20 ± 5 


To (AU) 


1 


T (K) (r=r D ) 


4500 ± 500 


q 


0.62 ± 0.03 


Disk mass (Mq) 


0.014 ± 0.001 


P 


1.3 ± 0.1 


Av tur b (km/s) 


0.8 ± 0.2 



Table 1 R Mon flat disk parameters. Note errors are those ob- 
tained from the fitting procedure and do not account for other 
uncertainties inherent of this kind of calculations (instrumental 
errors, optical depth effects, etc.). 



flat disk model. Our model assumes Local Thermo- 
dynamic Equilibrium (LTE) and standard 12 CO abun- 
dance equal to 8T0 -5 . With these assumptions, both, 
the CO 2— >1 and 1— >0 emissions are well fitted with a 
disk of 0.014 Mq in Keplerian rotation around the star. 
The disk radius is 1500 AU and the inclination angle is 
20° from an edge-on view. The mass of R Mon is fitted 
to 8 Mq . The gas kinetic temperature and density varies 
as T (K) = 4500 r"« and p (cm' 3 ) = Po r~P, with 
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Fig. 2 This figure shows the integrated intensity map of the 
12 CO 2^1 emission at blue (3.6-7 km/s) and red (12-15.4 km/s) 
velocities separately. The emission presents a bipolar morphology 
with a blue parabolic lobe towards the North and the red one 
towards the South. This N-S distribution is consistent with the 
morphology of the bipolar outflow at larger scales. The red lobe 
shows a jet-like morphology, i.e., we have detected a molecular 
jet associated with this BO star. The strip-disk line shows the 
expected orientation of the disk. 



q=0.62 and p=1.3. The density profile obtained with 
the model is consistent with the mass derived from the 
continuum data within a factor 2. A turbulent width 
component in the velocity profile is set to 0.8 km/s and 
added to the thermal width (table [1]) . 

In figure [4] the individual fit to some of the spectra 
(corresponding to certain position offset in AU along 
the major axis of the disk) for 12 CO 2— »1 is shown. A 
very good fit is obtained all over the disk, even in other 
directions outside the major axis. This is the strongest 
evidence of a disk in Keplerian rotation obtained so far 
around an early Be star. 

Since the 12 CO lines are optically thick, only the 
surface of the disk is contributing to the emission, so 
the mass of the gaseous disk cannot be derived from 
these observations. As a consequence, the assumed ge- 
ometry for the disk determines its mass. In fact, we 
reasonably fit the 12 CO emission using both a flat and 
a flared disk model, but with a mass 8 times larger in 
the second case. We need to observe the rarer isotope 
13 CO to determine the opacity of the CO lines and de- 
rive the mass and density structure of the disk, and 
eventually the disk geometry (see discussion in Paper 
I). In figure [5] we present the model results for 13 CO, 
showing that the flat disk model also fits this emission 
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Fig. 3 At the right the PV diagrams of the CO rotational lines 
along the strip drawn in Fig. [2] are shown. White contours super- 
imposed show the synthesized P-V diagrams obtained with out 
flat disk model. Contour levels are 200-800 mjy/beam in steps of 
100 mjy/beam (3.16 K) for the CO l->0 diagram, and 200-800 
mjy/beam in steps of 200 mjy/beam (3.46 K) for the CO 2— »1 
one. At the left, a comparison between the observed (grey) and 
the modeled (dark) spectra is shown for the star position. Note 
that the emission at the velocities of the molecular cloud (9.5 ± 
2.5 km/s) has been filtered out by the interferometric observa- 
tions. 
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Fig. 5 The 13 CO disk emission detected towards R Mon is well 
fitted by the flat disk model. The flared disk model (with a mass of 
0.08 Mq) predicts an emission clearly above the observed values. 



using the same parameters and an abundance of 9T0~ 7 
cm~ 2 (90 times below 12 CO standard abundance, as ex- 
pected). The flared disk model predicts an emission that 
exceeds the observed one by a factor >2. As a result of 
these new observations and modellings, now we can con- 
firm that the disk is flat. Furthermore, the 13 CO/ 12 CO 
line intensity ratio reveals as a power diagnostic tool to 
discern between flared and flat disks. 
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Fig. 4 These figures show the observed and the flat disk model predicted spectra over different positions along the major axis of the 
disk, for 12 CO J = 2-*l. Offsets in AU are labelled at the top right corner in each figure. The horizontal and vertical lines represents 
respectively the assumed noise level (3 K) and interferometric filtering limits (v ± 2.5 km/s,) in the fitting process. 

4 Conclusion References 



Based on IR observations several authors have proposed 
that, contrary to HAe and TT stars, HBe stars seem 
to have geometrically flat disks \Acke et al. 20 05). Our 
result in R Mon confirms this point. Moreover, the slope 
parameters for density and temperature distributions 
along the R Mon disk are similar to those found on TTs 
and some HAe \Dutrey et al. 20 04 )i an d, as most of 
them, is also in Keplerian rotation. The main difference 
respect to TTs and HAe stars is the mass of the disk, 
which is an order of magnitude lower in R Mon and in 
other HBe stars \Fuente et al. 2003\ . 

The similarities could indicate a common mecha- 
nism for star formation in the range 1—8 M Q . On 
the other hand, the different mass of the disks clearly 
indicates a shorter timescale for the evolution of disks 
around HBe stars, where the material is quickly dis- 
persed. In fact, the low occurrence of disks found around 
HBe stars \Fuente et al. 200~3\ supports this conclu- 
sion. A rapid evolution is also supported by the value 
of the spectral index (3. This indicates the presence of 
large dust grains, which drops the optical depth, so the 
UV radiation penetrates to photoevaporate the exter- 
nal layers (Dullemond and Dominik 2004 I- Such disks 



could evolve to flat geometries in 10 5 years. 

Future investigations will contribute to clarify this 
scenario. In this way, ALMA will be the definite instru- 
ment to study the properties of disks around interme- 
diate mass protostars. 
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